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We describe computer algorithms that produce the complete set of isohcdral tilings by 
n-omino or n-iamond tiles in which the tiles are fundamental domains and the tilings 
have 3-, 4-, or 6-fold rotational symmetry. The symmetry groups of such tilings are of 
types p3, p31m, p4, p4g, and p6. There are no isohedral tilings with p3ml, p4m, or 
p6m symmetry groups that have polyominoes or polyiamonds as fundamental domains. 
We display the algorithms' ou tput and give enumeration tables for small values of n. 
This expands earlier worksl^lAI 

Keywords: polyomino; polyiamond; isohedral tiling. 
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1. Introduction 



Polyominoes and polyiamonds and their tiling properties have been the subject of 
computational geometry research that investigated which polyominoes can tile the 
plane isohedrally and which can tile by translations alone. LU^. In earlier papers"^, 
we gave algorithms to create polyominoes and polyiamonds that were fundamental 
domains for isohedral tilings having p3, p4, or p6 symmetry groups. In this paper, 
we consider the expanded task of producing polyomino and polyiamond tiles that 
generate isohedral tilings of types p3, p3ml, p31m, p4, p4m, p4g, p6 or p6m 
and for which the tiles are fundamental domains of the tiling. Recently an extension 
of this study was published as a separate paper^ in which we carried out these 
investigations for symmetry groups of types pmm, pmg, pgg, and cmm. 

A polyomino (or n-omino) is a tile homeomorphic to a disk, made up of n 
unit squares that are connected at their edges; that is, the intersection of two unit 
squares in the polyomino is either empty or an edge of both squares. Similarly a 
polyiamond (or n-iamond) is a tile homeomorphic to a disk, made up of n unit 
equilateral triangles that are connected at their edges; the intersection of two unit 
triangles in the polyiamond is either empty or an edge of both triangles. 

An isohedral tiling of the plane is a tiling by congruent tiles in which the sym- 
metry group of the tiling acts transitively on the tiles. A fundamental domain for 
an isohedral tiling is a region of least area that generates the whole tiling when 
acted on by the symmetry group of the tiling. A fundamental domain for an iso- 
hedral tiling cannot contain two points that are identical under the action of the 
symmetry group of the tiling. This leads immediately to the following observations: 

(i) In an isohedral tiling in which each tile is a fundamental domain, no tile can 
contain a rotation center or axis of reflection for the whole tiling except on its 
boundary. 

(ii) In an isohedral tiling having a polyomino as fundamental domain, a J^-fold center 
of the tiling can only occur at a vertex of a square that is a "corner" of the 
polyomino; that is, only one unit square of the polyomino contains that vertex. 

(iii) In an isohedral tiling having a polyiamond as fundamental domain, 3-fold or 
6-fold centers of the tiling can only occur at vertices of unit triangles on the 
boundary of the polyiamond. Moreover, a 3-fold center can only occur at a vertex 
at which exactly two unit triangles meet. A 6-fold vertex can only occur at a 
vertex that is a "corner" of the tile, that is, exactly one unit triangle of the 
polyiamond contains that vertex. 

(iv) There are no 3-or 6-fold centers in a tiling by polyominoes. Hence, there are no 
p3, p3ml, p31m, p6 or p6m isohedral tilings by polyominoes. There are no 
4-fold centers in a tiling by polyiamonds. Hence, there are no p4, pA m , or VAQ 
isohedral tilings by polyiamonds. 

In each of the sections that follow, we begin with a fixed lattice of symmetry 
elements for a symmetry group G (that is, a fixed array of rotation centers, reflection 
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axes and glide-reflection axes for elements of G) and give a backtracking procedure 
to produce a complete set of polyominoes (or polyiamonds) that are fundamental 
domains for G. A tile T is a fundamental domain for a symmetry group G if the 
action of G on T produces an isohedral tiling and T is a region of minimal area 
for which G can generate that tiling. G will be contained in (or equal to) the full 
symmetry group of the tiling. 

As we consider symmetries of our isohedral tilings, the following theorem will 
be useful. "J 



Theorem 1. Let G be one of the 11 two-dimensional symmetry groups and T an 
isohedral tiling generated by G acting on a tile T that is a fundamental domain for 
G. Let G' be the full symmetry group of T ■ If G is a proper subgroup of G' , there 
is element of G' (other than the identity) that leaves T invariant. In this case, a 
fundamental domain for G' has area smaller than T. 



2. p4 

2.1. Creating polyominoes as fundamental domains for p4 
symmetry groups 

We begin with a lattice of unit squares as shown in Fig. [T] Our n-omino tiles that 
are fundamental domains for a p4 group will be built from these unit squares. By 
observation (ii) in section [TJ the 4-fold rotation centers for a p4 symmetry group 
must be located at lattice points. So, first we place a 4-fold rotation center, (a black 
circle) at a lattice point and call this the origin. We then place orthogonal unit 
vectors u and v at the origin (see Fig. [T]). Next we place a second 4- fold rotation 
center, the white circle in Fig.[TJ at xu + yv, where x and y are nonnegative integers, 
not both 0. These two choices of 4-fold rotation centers determine a whole p4 lattice 
of rotation centers; rotations about the two placed black and white centers generate 
a p4 group G. 





J~ 


f\ 


I i 

V 


II p 

u 



Fig. 1. A lattice of unit squares with black and white 4- fold rotation centers that generate a 
p4 group. The shaded isosceles right triangle is a standard fundamental domain; its unmarked 
third vertex is a 4- fold center equivalent to the white center. Here x = 1, y = 3; the area of the 
fundamental domain is 5 sq. units. 
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The area S of a fundamental domain for this p4 group G is given by 

s= x ^ (1) 

taking the area of a unit square as 1. Since we want our n-omino to be a fundamental 
domain, n = S and S must be an integer. Therefore, x and y must have the same 
parity, and 

n = 1, 2, 4, 5, 8, 9, 10, 13, 16, 17, 18, 20, 25, 25, 26, 29, 32, 34, 36, 37, 40, 41, 45, 49, . . . 

(2) 
where two pairs (x,y), namely, (5,5) and (7,1) correspond to n = 25. All unit 
squares in the lattice are classified into n equivalence classes by the action of the 
p4 group G. We denote the equivalence class of a unit square e as C(e). 

We construct a set 3F n of n-ominoes that are fundamental domains for the p4 
group G and a given n by following Procedure 1 below, using these definitions: 

(1) T is a set of unit squares; B(T) is a set of unit squares that are edge-adjacent 
to the squares in T; U n is a set of n-ominoes. 

(2) When T is the empty set, we define B{%) as the set of four squares around the 
origin (the four unit squares at the lower left in Fig. [TJ) 

(3) We define E(T), the boolean function of T, which is true if #T = n and the 
white circle is on the boundary of T. Otherwise E(T) is false. 

(4) B'{T) = {e|e € B{T),C{e) ^ C(/),for all/ e T}. This is the set of all unit 
squares that are edge-adjacent to those in T, but not equivalent to any unit 
squares in T. 

Procedure 1 creates a sequence of pairs of sets (T, Ut), in which Ut is the set of 
unit squares that can be added to T to create the next set T in the sequence. The 
sets T are built up, one element at a time, until a set T is achieved for which E(T) 
is true, at which time T is added to the collection S? n . When this takes place, the 
procedure backtracks to the most recent previous pair (T, Ut) for which Ut ^ 0, 
and repeats the process. 

Procedure 1. 

(1) Fix n (from the list in (2) above). Begin with 5^ empty. 

(2) Make T empty. Make U T = B'(T) = B{%). Make k = 0. Make S k = {(T, U T )}. 

(3) For (T, Ut) in Sk, if Ut ^ 0, choose an element e (a unit square) of Ut- Remove 
e from Ut and save the pair (T, Ut) in Sk- 

(4) Increase k by 1. Add e (from step 3) to T to create a new fc + 1-omino T, and 
let C/ T = B'(T) for the new T. Add the new pair (T, U T ) to 5 fc . 

(5) If T = and C/ T = 0, quit the procedure. 

(6) If E(T) is true, an n-omino tile is completed and add T to S? n provided there is 
no equivalent tile in 2? n . We regard two n-ominoes equivalent if conditions (a) 
and (b) below are satisfied. 



July 4, 2011 5:25 WSPC/Guidelines ws-ijcga-fkd3x 



Polyominoes and Polyiamonds as Fundamental Domains of Isohedral Tilings 5 

(a) The tiles are congruent (including mirror reflection). 

(b) When the tiles are superimposed, the positions of m-fold rotational centers 
on the boundaries of the tiles are the same. If there are several m-fold ro- 
tation centers for the same m, they can be permuted appropriately before 
comparison. 

(7) If U T = 0, decrease k by 1. 

(8) Go back to step 3. 

Fig. [2] shows the set of inequivalent n-ominoes in 2T n for n < 8. From each of 
these n-ominoes we can obtain the associated p4 tiling by using the black and white 
circles as 4-fold rotation centers. Fig. [3] shows the corresponding isohedral tilings 
produced by n-ominoes in Fig. [2] for n < 5. 



o^<V^>\\^>^. • 



(1-1) (2-1) 4-1 4-2 (4-3) 5-1 



5-1-2 (5-2) (5-2-2) 5-2-3 (5-3) 5-4 




(8-4) 



(8-6) 



18-7) 



(8-< )> 




Fig. 2. List of n-ominoes produced by the procedure in section 2.1 for n < 8. The tiles are 
fundamental domains for the p4 group used to construct them. The labels indicate n followed 
by the tile number for that n. Parentheses indicate that the tiles produce tilings having more 
symmetries than the p4 group that generates the tilings. 



2.2. Symmetries of tiles 

The list of isohedral tilings corresponding to tiles in Fig. [5] includes some tilings 
having symmetry group larger than the p4 group G generated by 4-fold rotations 
about the black and white rotation centers. By Theorem[TJ in every case where this 
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5-8 

Fig. 3. List of p4 isohedral tilings by n-ominoes in Fig. [2] for n < 5, generated by a given p4 
group. Labels correspond to those in Fig. [2] The symmetry group of each tiling is the given p4 
group, except for tilings whose labels are in parentheses. 



occurs, the n-omino that generates the tiling must have reflection and/or rotation 
symmetry. 

We outline below how to identify polyominoes whose corresponding isohedral 
tilings have symmetry groups larger than the group G that generated the tiling. 
Such tiles in Fig. [5] and their tilings in Fig. [3] are identified by parentheses around 
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their labels. 

• Select a polyomino that has rotation and/or reflection symmetry, and examine 
its tiling T generated by G. 

• Look at all vertices and centers of unit squares in a polyomino in T except for the 
original 4- fold centers we have chosen (black and white centers), and determine 
whether or not T is invariant under a 4-fold rotation about such a point. If so, 
then it is a new 4-fold center for 7", and this symmetry is not in G. For example, 
in Fig. |3j tilings (5-2-2) and (5-2-3), while generated differently, are the same 
tiling and have 4- fold centers at every vertex where 4 tiles meet. Tiling (5-3) has 
4- fold rotation centers at the centers of the "cross" tiles and at every vertex where 
4 tiles meet. These tilings have full symmetry groups of type p4, with G as a 
proper subgroup. 

• Other new symmetry elements of T can be sought by using chart 2 irM If the line 
joining the black and white 4-fold centers we placed is an axis of mirror reflection 
for T, then T is type p4m. For example, in Fig. [3J tilings (1-1) and (4-3) are 
type p4m. If this line is not a mirror reflection axis, but the line connecting two 
adjacent (nearest) 2-fold centers for T is a mirror reflection axis for T, then 7~~ is 
type p4g. For example, in Fig. [31 tiling (2-1) is type p4g. 

In sections [3] and 2] that follow, we indicate how to modify our Procedure 1 so 
as to produce directly polyominoes that are fundamental domains for p4g or p4m 
symmetry groups. 

3. p4g 

3.1. Creating polyominoes as fundamental domains for p4g 
symmetry groups 

A p4g symmetry group contains 4-fold rotations, reflections, and glide-reflections; 
the subgroup generated by its 4-fold rotations is type p4. Thus we begin as in 
section [21 with a lattice of unit squares. Following observation (ii) in section [1] we 
place a 4-fold rotation center (a black circle) at a lattice point and call this the 
origin. Orthogonal unit vectors u and v are then placed at the origin (see Fig. [4j. 

Next, we place the reflection axes that are nearest to the origin; these must lie 
along the edges of unit squares according to observation (i) in section [TJ These axes 
are vertical and horizontal lines at iu and at iv, respectively, where X is a positive 
integer. The placement of the origin and choice of a: determine the whole p4g lattice 
of rotation centers, reflection axes and glide-reflection axes since the p4g group G is 
generated by 4- fold rotations about the origin and reflections in one of the reflection 
axes nearest the origin. 

The area S of a fundamental domain for the group G is given by 

S = x 2 (3) 

where the area of a unit square is 1. Since we want our n-omino to be a fundamental 
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ifc; 



s-'— 6-- 



<?- 



Fig. 4. A lattice of unit squares with symmetry elements for a p4g group. Black and white circles 
are 4- fold rotation centers, thick solid lines are reflection axes, and dashed lines are glide-reflection 
axes. 2-fold rotation centers are at the intersections of reflection axes. 4-fold rotations about the 
center black circle (the origin) and reflections about one axis nearest the origin generate the p4g 
group. A standard fundamental domain is shaded. Here x = 2; the area of a fundamental domain 
for the p4g group is 4. 



domain, n = S. Therefore, 



1,4,9,16,25,36,49,64, 



(4) 



Action by the group G partitions the unit squares into n equivalence classes; as 
before, we denote the equivalence class of a unit square e as C(e). We construct 
the set ST n of n-ominoes that are fundamental domains for G using Procedure 1 in 
section [2l with n chosen from the list in ([4]) and the additional constraint that unit 
squares in the set T must be in the region bounded by the mirror reflection axes 
nearest to the origin. 

Fig. [5] shows the set of incquivalent tiles in !j n for n < 9. We can obtain the 
associated p4g tilings by using the origin as a 4-fold rotation center to fill out 
the square bounded by the reflection axes nearest the origin, then reflecting this 
square in its edges. Fig. [6] shows the corresponding isohedral tilings produced by 
polyominoes in Fig. [5l 



□ 



(1-1)4-1 (4-2) 4-3 9-1 









9-14 (9-15) 9-16 9-17 9-18 9-19 9-20 9-21 9-22 



Fig. 5. List of n-ominoes produced as described in section 3.1 for n < 9; these are fundamental 
domains for the p4g group used to construct them. The labels indicate n followed by the tile 
number for that n. Parentheses indicate that the tiles produce tilings having more symmetries 
than the p4g group that generates the tilings. 
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Fig. 6. List of p4g isohedral tilings by n-ominoes in Fig. \E\ generated by a given p4g group. 
Labels correspond to those in Fig. [5] The symmetry group of each tiling is the given p4g group, 
except for tilings whose labels are in parentheses. 
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9-25 
Fig. W\ (Continued) 



9-26 



3.2. Symmetries of tiles 

The only tilings in Fig. [6] having symmetry group larger than the p4g group G that 
generated them are tilings with square polyomino tiles. In fact, this is always true. 



Theorem 2. Every tiling by a polyomino T produced by our algorithm in section \37l\ 
has as its full symmetry group the p4g group that generated it and T is a funda- 
mental domain for the tiling, except in the case when the shape of T is square. In 
that case, the tiling has symmetry group p^m and T is not a fundamental domain. 



Proof. Let G' be the full symmetry group of a tiling T produced by a p4g group 
G acting on a polyomino T that is a fundamental domain for G, as described in 
section IBTTl If G is a proper subgroup of G", then by Theorem [TJ a rotation or 
reflection symmetry in G' is also a symmetry of T ■ By its construction, one corner 
of T is a 4-fold center for G at the origin; without loss of generality, we may assume 
that T contains the unit square shown in black in Fig. [7] T can also be obtained 



■==i zzz^wjzz zzz_mjzz ---zm 



(a) 



(b) 



(c) 



(d) 



Fig. 7. The derivation of polyomino fundamental domains T for a p4g group by removing unit 
squares from and adding equivalent unit squares to the original square polyomino fundamental 
domain in (a). Here x = 4. 
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beginning with anixi square polyomino tile as in Fig. [TJa), then removing unit 
squares and adding other unit squares equivalent (by a rotation of 90° or 270° 
about the origin) to these, always keeping the new tile homeomorphic to a disc. 
Figures EJb), (c), and (d) illustrate this. From this, we see that one edge of T that 
abuts a reflection axis in G has at least x unit squares in a straight row. Since the 
two edges of T that join the origin to each of the reflection axes at a;u and xv are 
related by a 90° rotation and are not straight unless T is a square, the only possible 
T that can have rotation symmetry is the square polyomino in Fig. [TJa). 

Suppose that T is not square but has reflection symmetry. If the two edges of 
T that abut the reflection axes at xu and xv are the same length, that reflection 
must map those edges to each other. But that diagonal reflection will not leave T 
invariant (see Fig. E]Jb)). If the two edges of T that abut the reflection axes at xu 
and a;v have different lengths, a reflection symmetry of T must leave fixed the longer 
edge, so must have its axis a horizontal or vertical line. But this reflection cannot 
be a symmetry in G' . For example, if X" is the image of T rotated by 90° about the 
origin, this reflection will not map T' onto another tile in T (see Fig. OJd)). □ 



4. p4m 

The symmetry elements of a p4m group are shown in Fig. [8j Black and white 



Fig. 8. The symmetry elements of a p4m group. Black and white circles are 4-fold centers, solid 
lines are reflection axes and dashed lines are glide-reflection axes. Points where glide-reflection 
axes intersect are 2-fold rotation centers. The shaded area is a fundamental domain for the p4m 
group. 



circles are 4-fold centers, solid lines are mirror reflection axes, and dashed lines are 
glide-reflection axes. Points where glide-reflection axes intersect are 2-fold centers. 
The shaded region is a fundamental domain for the p4m group. A polyomino that 
is a fundamental domain for the tiling must have its edges on the reflection axes, 
by observation (i) in section Q] But this is clearly impossible. 

Theorem 3. There are no p^m isohedral tilings having polyominoes as fundamen- 
tal domains. 
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5. p3 

5.1. Creating polyiamonds as fundamental domains for p3 
symmetry groups 

To build our n-iamond tiles that will be fundamental domains for a p3 isohedral 
tiling, we begin with a lattice of unit equilateral triangles. By observation (iii) in 
section [TJ the 3-fold rotation centers for a p3 symmetry group are located at the 
lattice points. So first we place a 3-fold rotation center, a black circle, at a lattice 
point and call this the origin, then place vectors u and v at the origin along edges 
of a unit triangle. Next we place a second 3-fold rotation center, a white circle, at 
xu + yv, where x and y are nonnegative integers, not both zero. See Fig.[9ta). These 




(a) (b) 

Fig. 9. (a) A lattice of unit triangles with two 3-fold rotation centers that generate a p3 group. 
Here X = 2, y = 1. (b) The shaded region is a standard fundamental domain for the p3 group 
generated by the black and white 3-fold centers; in this example, the area of a fundamental domain 
is 14 triangular units. 



two choices of 3-fold rotation centers determine the whole p3 lattice of rotation 
centers; 3-fold rotations about the black and white centers generate the whole p3 
symmetry group. 

A standard fundamental domain for the p3 group G generated by the black 
and white 3-fold centers is a 120° rhombus having as short diagonal the segment 
that joins the black and white circles, as shown in Fig. HJb). The other two vertices 
of this rhombus (gray in Fig. Mb)) are also 3-fold centers for G. The area S of a 
fundamental domain for G is given by 

S = 2(x 2 +y 2 + xy) (5) 

taking the area of a unit triangle as 1. Since we want our n-iamond to be a funda- 
mental domain, n — S. Therefore, 

n = 2, 6, 8, 14, 18, 24, 26, 32, 38, 42, 50, 54, 56, 62, 72, 74, 78, 86, 96, 98, 98, ... , (6) 

where two pairs (x, y), namely, (5, 3) and (7, 0) correspond to n = 98. 

All unit triangles are classified into n equivalence classes by the action of G, and 
we denote the equivalence class of a unit triangle e as C(e). 
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We construct a set 2F n of n-iamonds that are fundamental domains for G by 
modifying Procedure 1 in section [21 choosing n from the list in ([6]), and replacing 
definitions in section [2] with these: 

(1) T is a set of unit triangles; B(T) is a set of unit triangles that are edge-adjacent 
to the triangles in T; 2F n is a set of n-iamonds. 

(2) When T is the empty set, we define B(%) as the set of six triangles around the 
origin in Fig.jQja). 

(3) We define E(T), the boolean function of T, which is true if #T = n and the 
white circle is on the boundary of T . Otherwise E(T) is false. 

(4) B'(T) = {e|e € B(T),C(e) ^ C(/),for all/ e T}. This is the set of all unit 
triangles that are edge-adjacent to those in T, but not equivalent to any unit 
triangles in T. 

The procedure creates a sequence of pairs of sets (T, Ut), in which Ut is the set of 
unit triangles that can be added to T to create the next set T in the sequence. 

Fig. [TO] shows the set of inequivalent n-iamonds in ^ for n < 8. From each of 
these we can obtain the associated p3 tiling by using the black and white circles as 
3- fold rotation centers. Fig. [11] shows the corresponding isohedral tilings produced 
by n-iamonds in Fig. [TO] 




(2-1) 



(6-1) 



(6-2) 



6-3 



(6-4) 




8-3 



8-4) 



8-5 



8-6 



8-7 



Fig. 10. List of n-iamonds produced as described in scction R-H for n < 8. The tiles are fundamental 
domains for the p3 group used to construct them. The labels indicate n followed by the tile number 
for that n. Parentheses indicate that the tiles produce tilings having more symmetries than the 
p3 group that generates the tilings. 



5.2. Symmetries of tiles 

The list of isohedral tilings in Fig. [TT] includes some tilings having symmetry group 
larger than the p3 group G generated by 3-fold rotations about the black and white 
rotation centers. By Theorem [TJ when this occurs, the n-iamond that generates the 
tiling must have reflection and/or 3-fold rotation symmetry. The full symmetry 
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(8-4) 



Fig. 11. List of p3 isohedral tilings by n-iamonds in Fig. 1101 generated by a given p3 group. 
Labels correspond to those in Fig. 1101 The symmetry group of each tiling is the p3 group that 
generated it, except for tilings whose labels are in parentheses. 



group of the tiling could be of any of these types: p3, p3ml, p31m, p6, and p6m. 

The tilings in Fig.[TTJthat have additional symmetries, and the n-iamonds in Fig.fTOl 
that generate them are indicated by parentheses in their labels. We outline below 
how to identify these polyiamonds. 



Select a polyiamond T that has rotation and/or reflection symmetry, and examine 
its tiling T generated by G. 

Look at all vertices and centers of unit triangles in T except for those centers 
of rotation in G (black, white and gray centers), and determine whether or not 
any can be new 3-fold centers for T. If a new 3-fold center is found, then the 
full symmetry group G' of T contains G as a proper subgroup and T is not a 
fundamental domain for T ■ If a new 3-fold center is a 6-fold center, G' is type 
p6m if there is a reflection axis joining it to any nearest 3-fold center; otherwise, 
G' is type p6. Tiling (6-4) in Fig. [IT] is t y P e p6m, and has a 6-fold center at the 
center of the hexagonal 6-iamond. If there are new 3-fold centers but none are 
6-fold centers, then G' is type p3. 

If the tiling contains only 3-fold centers for G, we seek new symmetry elements 
by using chart 2 irM 
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(a) If some 3-fold center for G is a 6-fold center, the tiling T will have p6m 
symmetry if the line joining a 6-fold center to a nearest 3-fold center is a 
reflection axis for T, otherwise it will have p6 symmetry. In Fig. Ill) tilings 
(2-1) and (8-4) have p6m symmetry, and tiling (6-1) has p6 symmetry. 

(b) Otherwise, if the line connecting two adjacent different 3-fold centers (say, 
black to white) is a reflection axis, T has p3ml symmetry. 

(c) Otherwise, if the line connecting adjacent 3-fold centers of the same kind 
(black-black, white-white, or gray-gray) is a reflection axis, T has p31m 
symmetry. In Fig. QTJ tiling (6-2) has p31m symmetry. 

6. p31m 

6.1. Creating polyiamonds as fundamental domains for p31m 
symmetry groups 



We begin with a lattice of unit triangles as in section 15.11 since a p31m symme- 
try group contains 3-fold rotations, reflections, and glide-reflections; the subgroup 
generated by its 3-fold rotations is type p3. To build an n-iamond tile that is a fun- 



Fig. 12. A lattice of unit triangles with a 3- fold rotation center at the origin (black circle), 
reflection axes (thick solid lines) and glide-reflection axes (dashed lines) that are symmetries in a 
p31m group. The 3-fold rotation about the black circle and reflections about one axis generate the 
p31m group, which also has 3-fold centers where reflection axes intersect (one of these is shown 
as a white circle). A typical fundamental domain for the group is shaded. Here x = 2; the area of 
the fundamental domain is 12 triangular units. 



damental domain for a p31m isohedral tiling, we first place a 3-fold rotation center, 
a black circle, at a lattice point and call this the origin. Then we place vectors u 
and v at the origin along edges of a unit triangle (see Fig. H"2"j) . 

Next, we place reflection axes that are nearest the origin; these must lie along 
the edges of unit triangles according to observation (i) in section [TJ These three 
axes intersect at the points x(u + v), x{— 2u + v), and x(u — 2v), where £ is a 
positive integer; their intersections are equivalent 3-fold rotation centers (one of 
these is shown as a white circle in Fig. [T2l . The placement of the origin and choice 
of x determine the whole p31m lattice of rotation centers, reflection axes and glide- 
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reflection axes. The p31m group G is generated by 3-fold rotations about the origin 
and reflections in one of the three reflection axes. 

The area S of a fundamental domain for the tiling is one-third the area enclosed 
by the three reflection axes closest to the origin, which is 

3x 2 



S 



(7) 



where the area of a unit triangle is 1. Since we want our n-iamond to be a funda- 
mental domain, n = S. Therefore, 



3, 12, 27, 48, 75, 108, 147, 192, 243, . 



(8) 



The action of G classifies all unit triangles into n equivalence classes and we 
denote the equivalence class of a unit triangle e as C(e). We construct a set ^ n of 
n-iamonds that are fundamental domains for G as described in section [5Tl choosing 
n from the list in (8) and having the additional constraint that the unit triangles 
in the set T must be in the region bounded by the mirror reflection axes in Fig. 1121 

Fig. IT31 shows the set of inequivalent n-iamonds in & n for n < 12. Fig. [j"4l shows 
the corresponding p31m tilings obtained by performing a 3-fold rotation about the 
origin (black circle) to fill out the triangle bounded by reflection axes (Fig. H2|) . 
then reflecting this triangle in its edges. 




12-11 12-12 12-13 12-14 12-15 12-16 12-17 12-18 12-19 



12-20 



Fig. 13. List of n-iamonds produced as described in section [6. II for n < 12. The tiles are funda- 
mental domains for the p31m group used to construct them. The labels indicate n followed by 
the tile number for that n. 



6.2. Symmetries of tiles 

Theorem 4. Every isohedral tiling by a polyiamond T produced by our algorithm 
in section \6.1\ has as its symmetry group the p31m group that generated it and T 
is a fundamental domain for the tiling. 



Proof. Our proof is analogous to that for Theorem [5J Let G' be the full symmetry 
group of a tiling T produced by a p31m group G acting on a polyiamond T that is 
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12-8 



12-9 



12-10 



12-16 



12-17 



12-18 




12-11 




12-19 



12-20 



Fig. 14. List of p31m isohedral tilings by n-iamonds in Fig. 1131 generated by a given p31m group, 
for n < 12. Labels correspond to those in Fig. 1131 Every ra-iamond in Fig. 1131 is a fundamental 
domain for its tiling in Fig. 1141 
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a fundamental domain for G, as described in section ["6. II If G is a proper subgroup 
of G' , then by Theorem [TJ a reflection symmetry or 3-fold rotation symmetry in G' 
is also a symmetry of T. T can be obtained from the shaded fundamental domain 
shown in Fig. [12] by removing unit triangles and adding unit triangles equivalent 
to these (by a 120° rotation about the origin), always keeping the new tile homeo- 
morphic to a disc. This process shows that one edge of T that lies on a reflection 
axis will have length at least x + 1 , and at most one other edge of T (lying on an 
adjacent reflection axis of G) will have this same length (see Fig. [j~3j) . From this, it 
follows that T cannot have 3-fold rotation symmetry. 

Suppose that T has reflection symmetry. If the two straight edges of T that lie 
on reflection axes for G have the same length, then the reflection that leaves T fixed 
must map these edges to each other. But this reflection cannot leave T invariant, 
since the other edges of T are related by a 3-fold rotation about the origin (see, for 
example, tiles 12-1, 12-2, 12-3 in Fig. IT3")l . If the two straight edges of T that lie on 
reflection axes for G have different lengths, then a reflection that fixes T must leave 
the longer edge fixed, and so is perpendicular to that edge. This can only happen 
if T is the rhombus fundamental domain in Fig. Q2J and in this case, the reflection 
symmetry for T is not a symmetry for the tiling T (see tiling 12-4 in Fig. |T4"|) . □ 

7. p3ml 

The lattice of reflections and 3-fold rotation centers of a p3ml symmetry group 
is shown in Fig. 115) black, white, and grey circles denote three inequivalent 3-fold 
rotation centers. Here, unlike the p31m case, all 3-fold centers lie on reflection 
axes. The shaded region bounded by reflection axes is a fundamental domain for 
the p3ml group that generates the tiling by reflections those axes. By observation 




Fig. 15. The lattice of reflection axes and 3-fold rotation centers of a p3ml group (glide-reflection 
axes are not shown). Three inequivalent 3-fold centers are black, grey, and white circles. Lines are 
reflection axes. The shaded region is a fundamental domain for the p3ml group generated by 
reflections in the axes surrounding the region. 



(i) in section [TJ this is the only polyiamond tile possible having the area of a fun- 
damental domain. But the full symmetry group of this tiling is type p6m, which 
has a fundamental domain with area 1/6 that of the shaded tile. Thus, 

Theorem 5. There are no p3ml isohedral tilings having polyiamonds as funda- 
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mental domains. 

We note that if the shaded region in Fig. [T3] is decorated with an asymmetric 
motif then there are n-iamonds (where n = k 2 , k a positive integer) having the 
shape of an equilateral triangle, for which the decorated triangle is a fundamental 
domain for a p3ml isohedral tiling. 

8. p6 

8.1. Creating polyiamonds as fundamental domains for p6 
symmetry groups 

A p6 symmetry group contains 6-fold, 3-fold, and 2-fold rotations; the subgroup 
generated by its 3-fold rotations is type p3. Thus we begin as in section EU with a 
lattice of equilateral triangles, to build n-iamond tiles that are fundamental domains 
for a p6 group. By observation (iii) in section[T] the 6-fold and 3-fold rotation centers 
for a p6 symmetry group are located at lattice points. So first we place a 6-fold 
rotation center, a black circle, at a lattice point and call this the origin, then place 
vectors u and v at the origin, along edges of a unit triangle. Next we place a 3-fold 
rotation center, a white circle, at xu + yv, where x and y are nonnegative integers, 
not both 0. See Fig.[16l These choices of 6-fold and 3-fold rotation centers determine 





(a) (b) 

Fig. 16. (a) A lattice of unit triangles with a black 6-fold rotation center and white 3-fold rotation 
center that generate a p6 group. Here x = 2, y = 1. (b). A standard fundamental domain for 
the p6 group generated by the black and white rotation centers; in this example, the area of the 
fundamental domain is 7 triangular units. 



the whole p6 lattice of rotation centers; rotations about the two chosen black and 
white centers generate the p6 group. 

The shaded triangular region shown in Fig. 116( b) is a standard fundamental 
domain for the p6 group generated by 6-fold rotations about the black circle and 
3-fold rotations about the white circle. Hence the area of a fundamental domain for 
this p6 group is given by 

S = x 2 + y 2 + xy (9) 

taking the area of a unit triangle as 1. Since we want our n-iamond to be a funda- 
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-\ 



/ 



A 



(1-1) 3-1 (4-1) 4-2 4-3 (7-1) 7-1-2 7-2 7-3 7-3-2 7-4 7-5-2 7-6 



7-7 7-8 7-9 7-10 7-11 7-12 7-12-2 7-13 7-14 7-15 7-16 9-1 



9-2 



9-3 (9-4J 9-5 9-6 9-6-2 9-7 9-8 9-9 



9-10 9-11 




9-11-2 9-12 9-13 9-14 9-15 9-16 9-17 



9-21 9-22 9-23 9-24 9-25 9-26 



Fig. 17. List of n-iamonds produced as described in section fe-H for n < 9. The tiles are fundamental 
domains for the p6 group used to construct them. The labels indicate n followed by the tile number 
for that n. Parentheses indicate that the tiles produce tilings having more symmetries than the 
p6 group that generates the tilings. 



mental domain, n — S. Therefore, 

n = 1, 3, 4, 7, 9, 12, 13, 16, 19, 21, 25, 27, 28, 31, 36, 37, 39, 43, 48, 49, 49, 



(10) 



where two pairs of (x, y), namely, (5, 3) and (7, 0) correspond to n = 49. 

The action of the p6 group classifies all unit triangles into n equivalence classes 
and we denote the equivalence class of a unit triangle e as C(e). We construct a 
set 3? n of n-iamonds that are fundamental domains for the given p6 group as in 
section l5Tl with n chosen from the list in (fTU|). Fig. [T71 shows the set of inequivalent 
n-iamonds in & n for n < 9. From each of these n-iamonds we obtain the associated 
p6 tiling by using the black circles as 6-fold rotation centers and white circles as 
3-fold rotation centers. Fig. [T8l shows the corresponding isohedral tilings produced 
by n-iamonds in Fig. 1171 



8.2. Symmetries of tiles 

Some tilings in Fig. [18] have symmetries in addition to those in the p6 group G that 
generated them. These tilings and their polyiamond tiles in Fig. [T7] are indicated 
by parentheses in their labels. We can identify these polyiamonds as follows. 

• Select a polyiamond that has rotation and/or reflection symmetry, and examine 
its tiling T generated by G. 
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Fig. 18. List of p6 isohedral tilings by n-iamonds in Fig. 1171 for n < 7, generated by a given p6 
group. Labels correspond to those in Fig. 1171 The symmetry group of each tiling is the given p6 
group, except for tilings whose labels are in parentheses. 



• If the line connecting adjacent 6-fold and 3-fold centers (black to white) is a 
reflection axis for T, the tiling has p6m symmetry. Tilings (1-1) and (4-1) in 
Fig. [TH] have p6m symmetry. These tilings also have 3- fold centers for the tiling 
at the centers of the polyiamond tiles. 

• Otherwise, look at all vertices and centers of unit triangles in a polyiamond in 
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7-12-2 



7-13 



7-14 



7-15 




7-16 
Fig. I J SI (Continued) 

T except for those black and white rotation centers in the p6 group we have 
generated and determine whether or not they can be new 3-fold centers for T. 
If new 3-fold centers are found, the full symmetry group of the tiling contains G 
as a proper subgroup. Tiling (7-1) in Fig. [18] has a larger p6 symmetry group 
than G since there are 3-fold rotation centers for T at the centers of the rotorlike 
polyiamond tiles. Note that the same tiles appear in tiling 7-1-2, but in this tiling, 
the centers of the tiles are not 3-fold centers of rotation for T ■ 

9. p6m 

Fig. [19] shows the symmetry elements of a p6m group. Black, white, and gray 



v '/f>x?>rv /^vOv 5 ' 



Fig. 19. The lattice of symmetry elements for p6m. Solid lines are reflection axes and dashed 
lines are glide reflections axes. Black and white circles are 6-fold and 3-fold centers, respectively, 
and smaller grey circles where two reflection axes intersect are 2-fold centers. The shaded region 
is a fundamental domain. 



circles are the 6-, 3-, and 2-fold centers, respectively. Solid lines are reflection axes 
and dashed lines are glidereflection axes. The shaded 30°- 60°- 90°-triangle is a 
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fundamental domain. A polyiamond that is a fundamental domain for the tiling 
must have its edges on the reflection axes, by observation (i) in section [T] But this 
is clearly impossible. Thus, 

Theorem 6. There are no p6m isohedral tilings having polyiamonds as fundamen- 
tal domains. 

10. Enumeration Tables 

Tables 1-3 give the number of tiles and isohedral tilings that we have generated in 
sections [2] through [9l N n is the number of inequivalent tiles T in 2F n for a given 
symmetry group G of type p3, p31m, p3ml, p4, p4m, p4g, p6, and p6m that 
generate an isohedral tiling having the n-omino or n-iamond tiles as fundamental 
domain. Tiles are equivalent only if they generate the same tiling by the action of 
the group G when each tile is marked with an asymmetric motif. For example, the 
tiles 5-2 and 5-2-2 in Fig. [3] are congruent and their corresponding tilings are the 
same, but the placement of their 4-fold centers is different, and so if the tiles are 
marked with an asymmetric motif, they generate different isohedral tilings. 

Marking each tile with an asymmetric motif also guarantees that the group 
that generates the tiling is the full symmetry group of the tiling. So N n is also the 
number of isohedral tilings having G as full symmetry group and having n-omino or 
n-iamond tiles as fundamental domain, when each tile is marked with asymmetric 
motif. 

S n is is the corresponding number of the N n tilings when the asymmetric motif 
of each tile is removed. That is, S n is the number of isohedral tilings having full 
symmetry group G and having (unmarked) n-ominoes or n-iamonds as fundamental 
domains. In our figures that depict the isohedral tilings for small values of n, these 
tilings do not have parentheses around their labels. These are the most important 
"counting" results in this article. N' n is the number of non-congruent tiles in the 
N n tilings, counted by ignoring rotation centers attached to tiles; similarly, S' n is 
the number of non-congruent tiles in the S n tilings. 



Table 1. p3, p31m and p3ml. 
p3 n-iamonds p31m n-iamonds p3ml n-iamonds 



n 


2 
1 


6 
4 


8 
7 


14 
306 


n 


3 


12 


n 


1 


4 


9 ... 


N n 


N n 


1 


20 


N n 


1 


1 


1 


J ,-, 





1 


6 


294 


O n 


1 


20 


»-> n 








... 


N' 

n 


1 


1 


7 


288 


N' 

n 


1 


20 


N' 
n 


1 


1 


1 


S" 





1 


6 


277 




1 


20 










... 



For example, for symmetry group p6 and n — 7, there are 20 tiles in Fig. 1171 
with corresponding isohedral tilings in Fig.[THl so N7 — 20. Since tiling 7-1 of Fig.H^l 
has parentheses around its label, Sj = 19. From Fig. [T71 we can see that four pairs 
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p4 n-ominocs 



Tabic 2. p4, p4g and p4m. 
p4g n-ominocs 



p4m n-ominocs 



n 


1 

1 


2 
1 


4 
3 


5 

12 


8 
15 


9 

82 


10 
300 


n 

N„ 


1 
1 


4 
3 


9 

26 


16 

596 


n 


1 


2 ... 


N n 


N n 





... 


*■ J n 








2 


9 


38 


77 


296 


O n 





2 


25 


595 


On 





... 


A" 
n 


1 


1 


3 


8 


15 


80 


277 


N' 

n 


1 


3 


26 


596 


N' 

n 





... 










2 


7 


38 


76 


275 


qi 





2 


25 


595 


"n 





... 



p6 n-iamonds 



Tabic 3. p6 and p6m. 

p6m n-iamonds 



n 



1 



7 



12 



N n 1 1 3 20 29 195 

S n 1 2 19 28 194 

N' n 1 1 3 16 27 191 

S' 1 2 16 26 190 



n 


1 


2 


3 ... 


N n 








... 


O n 








... 


N' 








... 


qi 








... 



of tiles are congruent: 7-1 and 7-1-2; 7-3 and 7-3-2; 7-5 and 7-5-2; 7-12 and 7-12-2. 
Thus Nij = 20 — 4 = 16. Among the 19 tiles counted for SV, there are also 16 
non-congruent tiles, so S' 7 = 16. 



11. Summary 

We have described computer algorithms that can enumerate and display isohcdral 
tilings by n-omino or n-iamond tiles for given n in which the tiles are fundamental 
domains and the tilings have 3-, 4-, or 6-fold rotational symmetry. Their symmetry 
groups are of types p3, p31m, p4, p4g, and p6. We have shown that there are 
no isohedral tilings with symmetry groups of types p3ml, p4m, or p6m that 
have polyominoes or polyiamonds as fundamental domains. For symmetry groups 
of types p3, p31m, p4, p4g, and p6 we used the backtracking Procedure 1 to obtain 
a set c3^ of n-omino or n-iamond tiles where each tile produced one isohedral tiling, 
generated by a given symmetry group G of one of these five types. We can denote 
^ n (G) as the set 2? n for that symmetry group G and ^*(G) the corresponding set 
of isohedral tilings. 

We investigated the symmetries of tilings in the set ^*(G) and noted those 
tilings that satisfy the following two conditions: (1) the full symmetry group of the 
tiling is G, and (2) the tiles are fundamental domains for G. We denote the subset 
of ,^*(G) that satisfies (1) and (2) as S^*{G) and the corresponding set of tiles as 
S^ n (G). (For small values of n, these tiles and their tilings were displayed with labels 
without parentheses.) The enumeration of J?* (G) is the main counting result of this 
article. Although the n-omino or n-iamond tiles produced by our algorithm are not 
always fundamental domains for the isohedral tilings they generate, if we mark these 
tiles with an asymmetric motif, then the set 3?* (G) is the set of all isohedral tilings 
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with symmetry group G in which the corresponding tiles in 2? n (G) are fundamental 
domains. The set 2?*(G) can then also include a marked fundamental domain for 
a p3ml symmetry group. In Tables 1-3 of section ITOl we used the notation N n — 
#^ n *(G)andS„ = #^„*(G9. 
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